Current lookback time-redshift bounds on dark energy 
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We investigate observational constraints on dark energy models from lookback time (LT) estimates 
of 32 old passive galaxies distributed over the redshift interval 0.11 < z < 1.84. To build up our 
LT sample we combine the age measurements for these 32 objects with estimates of the total age of 
the Universe, as obtained from current CMB data. We show that LT data may provide bounds on 
the cosmological parameters with accuracy competitive with type la Supernova methods. In order 
to break possible degeneracies between models parameters, we also discuss the bounds when our 
lookback time versus redshift sample is combined with with the recent measurement of the baryonic 
acoustic oscillation peak and the derived age of the Universe from current CMB measurements. 

PACS numbers: 98.80.-k; 95.36.+X 



I. INTRODUCTION 

The ratio of the dark energy pressure to its energy 
density, the so-called equation of state (EoS) parame- 
ter, to, is nowadays one of the most searched numbers 
in general relativistic cosmology (see, e.g., Ref. [l[ for 
recent reviews). This is so for at least two diferent rea- 
sons. First, because if one could set Li to be exactly — 1, 
then there would be a great probability of identifying the 
dark energy with the vacuum state of all existing fields 
in the Universe, i.e., the cosmological constant (A). Sim- 
ilarly, if a value u> =/= — 1 (or a time variation of to over 
the cosmic evolution) is unambiguously found, then one 
could not only rule out the cosmological constant but 
also seriously think of the dark pressure responsible for 
the current cosmic acceleration as the potential energy 
density associated with a dynamical scalar field <p. The 
possibility to ^ —I still leads to two diferent routes, i.e., 
a quintessence field if —1 < to < 1/3 [|J or a phantom 
component for lu < — 1 [3j . Both cases violate the strong 
energy condition, p + 3p > 0, but the latter goes further 
and also violates the null energy condition, i.e., p + p > 

1- 

Clearly, constraining the value to from diferent sets of 
observational data constitutes an important way to im- 
prove our understanding of the actual nature of the dark 
energy. In this regard, finding new methods or reviv- 
ing old ones that could directly or indirectly quantify the 
amount of dark energy present in the Universe, as well 
as determine its EoS parameter, are important tasks for 
both theoretical and observational cosmologists. Bounds 
on to have been obtained from observations based on com- 
pletely distinct physics @, i, 0, 0, & GjJ [ID, 03, El (see 
also [l| and Refs. therein). This diversity of technics, as 
well as combinations among them, is particularly impor- 



tant to a more realiable determination of to, since diferent 
methods may constrain different regions of the parame- 
ter space and, therefore, may be complementary to each 
other. 

In this paper, by following the methodology presented 
m Refs. 14, 151 , we are particularly interested in deriving 
the current lookback time (LT) versus redshift bounds 
on the dark energy EoS and its density parameter from a 
sample of 32 passively evolving galaxies, as recently stud- 
ied by Simon et al. The absolute age for these 32 ob- 
jects was determined by fitting stellar population models 
and the sample includes observations from the Gemini 
Deep Deep Survey (GDDS) [l7} and archival data 0]. 
The same data, along with other age estimates of high- 
z objects, were recently used to reconstruct the shape 
and redshift evolution of the dark energy potential [Ty] , 
to place bounds on holography-inspired dark energy sce- 
narios 11811 . as well as to constrain models of modified 
gravity [191 ]. Here, we focus our analysis on two dark en- 
ergy models, namely, the standard ACDM scenario and 
a flat universe driven by non-relativistic matter (bary- 
onic + dark) and a smooth negative-pressure dark en- 
ergy component (wCDM). In order to better constrain 
the parametric spaces for these scenarios, we also com- 
bine LT data with the recent measurement of the bary- 
onic acoustic oscillation peak [ll| and the derived age of 
the Universe from current CMB measurements [Si . 



II. LOOKBACK TIME-REDSHFT TEST 
A. Theory 

The lookback time-redshift relation, defined as the dif- 
ference between the present age of the Universe (to) and 
its age (t z ) when a particular light ray at redshift z was 
emitted, can be written as 
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FIG. 1: The age-redshift data points. Left: Original data from Ref. [lfj. This sample corresponds to 32 old passive galaxies 
distributed over the redshift interval 0.11 < z < 1.84 and includes observations from the Gemini Deep Deep Survey (GDDS) [TtJ 
and archival data 0. Right: The LT sample. To obtain this sample we have combined the age measurements for these 32 
objects with estimates of the total age of the Universe, to bs = 13.7 ± 0.2 Gyr, as obtained from current CMB data @]. 



where H 1 = 9.78h 1 Gyr and h ranges in the HST key 



project la interval 0.64 < h < 0.8 
expression, 



[20| | . In the above 



H(p) = a/e(p) with E(p) = n m &- 3 + n x ar 3 ( 1+ul) (2) 

where H(p) = H(p)/H , the complete set of parameters 
is p = Qj, uj (j = m, A and x stand for matter, cosmolog- 
ical constant (lo = —1) and dark energy (lj < 0) density 
parameters, respectively), and the subscript denotes 
present-day quantities. 

To proceed further, let us now consider an object at 
redshift m whose the age t(zi) is defined as the difference 
between the age of the Universe at Zi and the one when 
the object was born zf, i.e., 
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or, equivalently, 



t(Zi) = t L (z F ) - t L (Zi). 



(3) 



(4) 



From the above expressions, we can define the observed 
lookback time to an object at Zi as [13] 



tT{zi;r) 



th{zF) ~ t(Zi) 
= [t°o S ~ t(Zi)] 
= t° " S ~ t(Zi) - 



n, (5) 



where ti stands for the so-colled delay factor, which ac- 
counts for our ignorance about the amount of time since 
the beginning of the structure formation in the Universe 
until the formation time (ij) of the object i [7[. 



We estimate the best-fit to the set of parameters p by 
defining the likelihood function 



C age oc exp [-xl ge {z;p,Ti)/2] , 



(6) 



where xlqe is g iven by 
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Here. 



= erf + cfobs , erf is the uncertainty in the indi- 
vidual lookback time to the i th galaxy of our sample and 
a t ob S stands for the uncertainty on the total expanding 

age of the Universe (t° hs ). The prior on the total age 
of the Universe is justified by the fact that quintessence 
scenarios that are able to explain age estimates of high-z 
objects may not be compatible with the total expand- 
ing age up to z — (and vice-versa) (see, e.g, [21| for a 
discussion) . 

Another important aspect concerns the delay factor 
t. Note that in principle there must be variations in 
the value of r for each object in the sample (galaxies 
form at different epochs). Differently from Ref. [l5j], in 
the present analysis the delay factor Tj for each object 
of the sample is assumed as a nuisance parameter, so 
that we marginalize over them. Note also that, although 
involving a n number of integrations, this marginalization 
may also be analytically obtained by defining a modified 
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FIG. 2: Left: Confidence contours at 68.3% and 95.4% in the parametric space f2 m — h from LT analysis. As physically 
expected, the larger the value of f2 m the smaller the value of the Hubble parameter that is allowed in order to fit the LT data. 
Middle: Likelihood function for fi m . From this calculation, we have found f2 m = 0.259 ± 0.030 at 95.4% (C.L.). Right: The 
Sl m — SI a plane arsing from our LT analysis. The best-fit favours a spatially open universe with Slk — 0.7. 



log-likclihood function x 2 , i.e., 

(/•oo roo 
J ...J dr^exp 



-oxL e (^p> T *) 



where 



A; 



V2a t 



^=t L { Zl - V )-[C s -t{Zi)\ 



and E is the second term of the rhs of Eq. ([7]). 



B. Data 



(8) 
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show, respectively, the original age estimates and trans- 
formed lookback time as a function of the redshift for the 
32 galaxies of Ref . [l6[ • 



III. RESULTS 

In this Section we discuss quantitatively how our LT 
sample may place bounds on the EoS and the density 
parameters describing the dark matter and dark energy 
through a statistical analysis of the data. We perform 
our analyses in the context of two different dark energy 
models, namely, the standard ACDM scenario (flat and 
general curvature) and a spatially flat universe driven by 
non-relativistic matter (baryonic + dark) and a negative- 
pressure dark energy component (wCDM), as described 
inEq. ©. 



In order to apply the method outlined above, we use 
age estimates of 32 old passive galaxies distributed over 
the redshift interval 0.11 < z < 1.84, as recently analized 
in Ref. [lfj. The total sample is composed by three sub- 
samples: field early- type galaxies from Ref. [12], whose 
ages were obtained by using SPEED models of Ref. [23| ; 
20 red galaxies from the publicly released Gemini Deep 
Survey (GDDS) P3 - Ref. HH re-analized the GDDS old 
sample by using a different stellar population models and 
obtained ages within 0.1 Gyr of the GDDS collaboration 
estimates - and the two radio galaxies LBDS 53W091 and 
LBDS 53W069 Q. 

To build up our LT sample, we combine the ages of the 
above galaxy sample with estimates of the total age of 
the Universe ig hs , according to Eq. (5). In our analysis, 
we assume ig hs = 13.7±0.2 Gyr, as obtained from a joint 
analysis involving current data of the most recent CMB 
experiments (WMAP, DASI, VSA, ACBAR, MAXIMA, 
CBI and BOOMERANG) @. In Panels (la) and (lb) we 



A. ACDM 

In Figures (2a)- (2c) we show the first results of our sta- 
tistical analyses. By fixing u> = — 1 in Eq. ([2]), Panel 2(a) 
shows contour plots (68.3% and 95.4% C.L.) in the fl m — h 
plane for the x1 ge given by Eqs. ©-([5]) plus a Gaussian 
prior on the Hubble parameter, h = 0.72 ± 0.08, as given 
by the final results of the HST key project [20]. As phys- 
ically expected, the larger the value of O m the smaller 
the value of the Hubble parameter that is allowed by the 
statistical analysis in order to fit the ti,{z) estimates. At 
95.4% (C.L.), we have found 0.18 < Q m < 0.23 or, equiv- 
alently, 0.74 < fl A < 0.77 (Q A = 1 - Q ro ). 

Panel 2(b) shows the likelihood function for the mat- 
ter density parameter. The dotted lines are cuts in the 
regions of 68.3% and 95.4% (C.L.). From this calcu- 
lation, we have found fl m = 0.259 ± 0.030 at 95.4% 
(C.L.), which is in good agreement with current fl m es- 
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FIG. 3: Results from joint analyses involving LT and BAO data. Left: The plane f2 m — h for a flat ACDM model. From this 
analysis we found fi m = 0.25 ±0.02 and h = 0.73±g;g§ at 2a level. Left: 68.3% and 95.4% confidence contours in the f2 m — Oa. 
The best-fit model clearly favours a spatially flat model with Q k ~ 0.07 and f2 m = 0.24 ± 0.02 and fi A = 0.69 ± 0.08 at 95.4% 
(C.L.). 



timates from CMB [6j and other independent [24[ re- 
sults. We also show in Panel 2(c) the f2 m — 51a space 
when the flat condition is relaxed, i.e., by adding the 
term f^a -2 = (1 — f2 in — fl\)a~ 2 to the E(p) function 
of Eq. ((2)). Differently to the results from current SNe 
la data (which prefer a spatially closed universe), the 
best- fit scenario is an open universe with ~ 0.7 and 
0.03 < O m < 0.1 and 0.15 < Q A < 0.32 at 95% (C.L.). 



1. Joint Analysis 

As well known, the acoustic peaks in the cosmic mi- 
crowave background (CMB) anisotropy power spectrum 
is an efficient way for determining cosmological parame- 
ters. Because the acoustic oscillations in the relativistic 
plasma of the early universe will also be imprinted on to 
the late-time power spectrum of the non-relativistic mat- 
ter [25|, the acoustic signatures in the large-scale clus- 
tering of galaxies yield additional tests for cosmology. In 
particular, the characteristic and reasonably sharp length 
scale measured at a wide range of redshifts provide an 
estimate of the distance-redshift relation, which is a geo- 
metric complement to the usual luminosity-distance from 
SNe la. Using a large spectroscopic sample of 46,748 lu- 
minous, red galaxies covering 3816 square degrees out to 
a redshift of z = 0.47 from the Sloan Digital Sky Suvey, 
Eisenstein et al. [ll[ have successfully found the peaks, 
described by the .4-parameter, i.e.[26j|. 



A 



1/2 



ZBAO 



ZBAO 



r 2 (z B Ao;p) 
E(zbao;p) 



1/3 



0.469( 



-0.35 



0.98 ' 



±0.017, 



(9) 



which can be used to constrain cosmological scenarios 
that do not have a large contribution of dark energy at 
early times [l2|. In the above expression, zbao = 0.35 
is the redshift at which the acoustic scale has been mea- 
sured, T(zbao; p) = / ZBA ° dz /H(zbao) is the dimension- 
less comoving distance to zbaOj an d we have taken the 
scalar spectral index n s — 0.95, as given in Ref. [f|. 

In Figure (3a) we show contour plots (68.3% and 95.4% 
C.L.) in the Q m — h plane for the LT ± BAO combination. 
Note that, relative to the results shown in Panel (2a), 
the allowed parameter space is now considerably reduced, 
with the 2ct bounds lying in the interval fi m = 0.25 ±0.02 



and h = 0.73 



+0.02 

-0.03- 



By allowing for arbitrary curvature 



and using a proper generalization of Eq. (9), Fig. (3b) 
shows the O m — f^A space for the joint LT ± BAO analysis. 
The allowed parameter space now is considerably reduced 
relative to the former case of Fig. (2c). In agreement with 
current CMB results [|| , the best-fit model for this joint 
LT + BAO analysis clearly favors a nearly flat universe 
with Q k ~ 0.07 and ft m = 0.24 ± 0.02 and tt A = 0.69 ± 
0.08 at 95.4% (C.L.). 



B. uCDM 

With the usual assumption that the effective EoS, 
w ~ J uj(z)Q, x (z)dz/Q x (z), is a good approximation for 
the wide class of dark energy scenarios, from now on we 
discuss the bounds from LT and LT±BAO data on uj. 

Figure (4a) shows the parametric space 51 m — co al- 
lowed at 68.3% and 95.4% (C.L.) from LT data only. 
Note that, although the matter density parameter is well 
constrained by these data, a large interval for w is still 
allowed. In particular, the best-fit model happens for 
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FIG. 4: Results for ojCDM model Left: The Q m — w space allowed by LT data. Left: LT + BAO bounds on the S7 m — uj plane. 
As discussed in Ref. fl5l ] , the tighter results found in this combined analysis reflect the complementarity between LT and BAO 
measurements. 



values of fi m ~ 0.1 and ui ~ —0.46. Note also that, al- 
though the physics behind LT and SNe la observations 
are quite different, LT constraints on the f2 m — uj plane 
are very similar to those obtained from SNe la measure- 
ments. This amounts to saying that combinations of LT 
+ SNe la data are not able to break possible degenera- 
cies on the plane On, — uj, and also that LT data may 
provide bounds on the cosmological parameters with ac- 
curacy competitive with SNe la methods (see Ref. [l5| 
for a discussion). 

A different result arises when the current BAO mea- 
surement at z = 0.35 is added to the analysis (Fig. 
4b). In this case, both uj and O m intervals are more 
tightly constrained, with the best-fit values given by 
O m ~ 0.27 and uj ~ -1.04. At 95.4% (C.L.) we also 
found 0.25 < O m < 0.29 and -1.21 < uj < -0.88. As 
discussed in Ref. [HI, the tighter results found in the 
combined analyses [Figs. (3b) and (4b)] reflect the com- 
plementarity between LT and BAO measurements, which 
in turn makes possible to break the degeneracies inherent 
to the parametric plane f2 ra — to. 



IV. FINAL REMARKS 

The recent accumulation of independent observational 
results has opened up a robust window for probing the 
behavior of the dark component responsible for the cur- 
rent cosmic acceleration. However, most of the methods 
employed to place limits on the dark energy EoS (uj) or, 
more generically, on the parametric space O m — uj, are 
essentially based on distance measurements to a particu- 
lar class of objects or physical rulers (e.g., SNe la, CMB, 
galaxy clusters, etc.). In this regard, it is also particularly 
important to obtain accurate and independent bounds on 
the physical behavior of the dark energy, as well as on the 



other main cosmological parameters, from physics relying 
on different kinds of observations. 

In this paper, by extending and updating the results 
of Ref. [Tg] , we have followed this direction and studied 
the current constraints on the parametric space Q m — uj 
from age measurements of high-z galaxies, as recently 
discussed in Ref. pjj]. By using a sample of 32 pas- 
sively evolving galaxies distributed over the redshift in- 
terval 0.11 < z < 1.84, we have transformed age mea- 
surements into LT estimates (by using the current values 
for the total age of the Universe from CMB data) and 
discussed quantitatively how these current age data may 
constrain the parametric spaces f2 m — and fi TO — uj. 
We have shown that LT data may provide bounds on 
the cosmological parameters with an accuracy compet- 
itive with SNe la methods [see, e.g., Figs. (3a) and 
(3b)]. Due to the complementarity between LT (age) 
and BAO (distance) measurements, our best results are 
obtained when joint analyses involving these two differ- 
ent observables are performed. By assuming uj = — 1 
(ACDM) and allowing for arbitrary curvature, we have 
found O m = 0.24 ± 0.02 and Q A = 0.69 ± 0.08 at 95.4% 
(C.L.), which clearly favours a nearly flat universe with 
f2k ~ 0.07. For a spatially flat model dominated by a 
negative-pressure component with a constant EoS uj, we 
have found 0.25 < O m < 0.29 and -1.21 < uj < -0.88 
(95.4% C.L.), which is again close to the so-called con- 
cordance scenario obtained from the usual distance-based 
combination of SNe + BAO + CMB data. 
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of Eq. (9) for ar- 
be written as A = 

2/3 

|f2 fc |r(z B Ao) ) , where 



A direct generalization 
bitrary curvatures can 

nil 2 



H(^AO) 1 ' 3 



z BAOvl n fcl 

the function J-{x) is defined as smh{x),x, and sin(a;) 
respectively, for open, flat and closed geometries. 



